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ABSTRACT 

Moderate resolution data for 40 quasi-stellar objects (QSOs) at z « 2 were combined with spectra of 
comparable resolution of 59 QSOs with redshifts greater than 1.7 found in the literature to form a large, 
homogeneous sample of moderate resolution (~1 A) QSO spectra. These spectra were presented and 
the statistics of the Lyman a forest were discussed in Paper I. In this analysis, we demonstrate that a 
proximity effect is present in the data, ie. there exists a significant (S.Str) deficit of lines at Zabs ~ Zem- 
Within 1.5 Mpc of the QSO emission redshift, the significance does depend on QSO luminosity, in 
accordance with the theory that this effect is caused by enhanced ionization of hydrogen in the vicinity 
of the QSO from UV photons from the QSO itself. The photoionization model of Bajtlik, Duncan, & 
Ostriker permits an estimate of the mean intensity of the extragalactic background radiation at the 
Lyman limit. We compare the results of this standard analysis with those obtained using a maximum 
likelihood technique. If the spectrum of the background is assumed to be identical to that of each 
individual QSO, and if this background is assumed to be constant over the redshift range 1.7 < z < 3.8, 
then the best fit value for J(i'o) is found to be 1.4^q g x 10~^^ ergs s~^ cm"^ Hz^^ sr~^, using QSO 
redshifts based on the Lya emission line. Systemic QSO redshifts based on the [OIII] A5007 emission 
line for 19 objects in our sample show an average redshift of ~400 km s~^ with respect to Lya emission. 
Using redshifts based on [OIII] or Mg II for the 35 objects for which they are measured and adding 
400 km s~^ to the remaining QSO Lya redshifts gives a lower value of J{vo), 7.OI4 4 x 10~^^ ergs 
s~^ cm~^ Hz^^ sr^^. This value is in reasonable agreement with the predictions of various models of 
the ionizing background based on the integrated QSO luminosity function. Allowing for the fact that 
individual QSOs have different spectral indicies which may also be different from that of the background, 
we use the standard methods to solve for the HI photoionization rate, F, and the parameters describing 
its evolution with redshift. The best fit value for the HI ionization rate we derive is 1.9ti o x 
s~^, in good agreement with models of the background which incorporate QSOs only. Finally, we use 
simulated Lyman a forest spectra including the proximity effect to investigate curve-of-growth effects 
in the photoionization model used in the analysis. We find that the presence of lines on the saturated 
part of the curve-of-growth could cause our estimates of the background intensity to be overestimated 
by a factor of two to three. This large absorption line sample and these techniques for measuring the 
background and understanding the systematics involved allow us to place what we believe are the firmest 
limits on the background at these redshifts. 

Subject headings: diffuse radiation — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

The spectra of quasi-stellar objects (QSOs) blueward 
of Lya emission show a prodigious number of absorption 
lines primarily due to Lya absorption by intervening neu- 
tral hydrogen along the line of sight to the QSO (Lynds, 
1971; Sargent et al. 1980; Weymann, Carswell, & Smith, 
1981). In models of these Lya systems, they are in pho- 
toionization equilibrium with a background ultraviolet ra- 
diation field. This radiation field is modeled as the inte- 
grated emission from QSOs and young galaxies (Bechtold 
et al. 1987; Miralda-Escude & Ostriker 1990; Madau 1992, 
Meiksin & Madau 1993, Madau & ShuU 1996, Haardt k 
Madau 1996, Fardal et al. 1998). 



The basic trend in the number density evolution of Lya 
absorbers (12 <^ log(N///) < 16 cm~^) at high redshift in 
a given QSO spectrum is a steep power law increase with 
redshift (Sargent et al. 1980). However, a decrease in the 
number of lines near the emission redshift of an individ- 
ual QSO relative to the number of lines expected from the 
power law distribution has been observed (Carswell et al. 
1982, Murdoch et al. 1986; Tytler 1987; Bajthk, Duncan, 
& Ostriker 1988, hereafter BDO; Lu, Wolfe, & Turnshek 
1991, hereafter LWT; Bechtold 1994, hereafter B94). The 
simplest explanation for this "inverse" effect is enhanced 
ionization of HI in the vicinity of the QSO by ultraviolet 
photons from the QSO itself. Thus, the name "proximity 
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effect" is also used. This interpretation, along with the as- 
sumptions about the spectrum of the background and the 
photoionization of the nearby intergalactic medium (IGM) 
by the QSOs, allows for a measurement of the mean in- 
tensity of the ionizing background at the Lyman limit of 
hydrogen (BDO), which can be compared to estimates of 
the integrated emission from QSOs. 

BDO find that J(i'o) is approximately constant over a 
redshift range 1.7 < z < 3.8. Expressing J(i'o) as J_2i 
X 10^^^ ergs s^^ cm^^ Hz^^ sr^^, these authors find 
J-21 ~ 1. Their best fit form for the dependence of J(i'o) 
on redshift rules out a luminosity dependent cutoff in the 
QSO luminosity function (BWLM; Schmidt, Schneider, & 
Gunn 1986). Several other authors have carried out this 
analysis on other data sets (LWT, Kulkarni & Fall 1993, 
hereafter KF93, B94, WiUiger et al. 1994, Cristiani et al, 
1995, Giallongo et al. 1996.) LWT found J_2i « 1 for 
1.7 < z < 3.8; and B94 found a value consistent with this, 
J-21 ~ 3 for 1.6 < z < 4.1. Srianand & Khare (1996) 
compile a sample of 69 QSOs from the literature (54 from 
B94) and obtain J_2i « 6 for 1.7 < z < 4.1. WiUiger et 
al. (1994) find a lower value, J_2i « 0.1 - 0.3 for z ~4.2 
from a single z=4.5 QSO; and Giallongo et al. (1996) find 
J_2i w 0.5 for 2.8 < z < 4.1, based on high resolution 
spectra of six QSOs. KF93 re-express the BDO formalism 
using maximum likelihood techniques and use this to de- 
rive J(fo) w 6 X lO^^'' ergs s~^ cm^^ llz~^ si~^ for z « 
0.5 from the sample of Bahcall et al. (1993). 

In this paper, a homogeneous sample of QSO spec- 
tra is used to measure the mean intensity of the ioniz- 
ing background via the standard proximity effect analysis 
and the maximum likelihood analysis of KF93. We have 
presented new Lyman a forest data for 39 objects with 
1.9 < Zem < 2.5 in Paper I; and supplemented this sample 
with 59 objects from the literature, and one QSO from 
the Hamburg/CfA Bright Quasar Survey (Dobrzycki, En- 
gels, & Hagen 1999). The spectra comprising our dataset 
are of moderate resolution, ^^1 A FWHM; but this is not 
a disadvantage, as this analysis requires good absorption 
line statistics and therefore many QSO sight lines. This 
is difficult to achieve at high resolution, the major reason 
for using a large set of moderate resolution spectra such 
as this one. 

A major uncertainty in the proximity effect analysis is 
in the systemic redshifts of the QSOs. Redshifts measured 
from low ionization permitted lines (e.g. Balmer lines or 
Mg II) or forbidden lines (e.g. [OIII] AA4959,5007) lines 
have been shown to be redshifted with respect to Lya and 
C IV emission by up to «250 km s^^ (Boroson & Green 
1992, Laor et al. 1995). B94 found that increasing the 
values of the QSO redshifts by 1000 km s~^ caused the 
best fit value of J_2i to be decreased by a factor of 3. We 
therefore obtained emission line spectra for several objects 
in our sample in order to examine redshift differences be- 
tween Lya and [OIII] AA4959,5007, Mg II or Balmer emis- 
sion. We investigate the effect of these shifts on the value 
of J(i'o) derived. 

In Section ^ the spectra used to measure Lyman limit 
fluxes for some of our sample objects and spectra used to 
perform systemic redshift measurements of several QSOs 
in our sample are presented. In Section || the techniques 
used to measure the mean intensity of the HI ionizing back- 



ground at the Lyman limit are discussed, which inclu des 
discussions of the QSO systemic redshifts in Section 
and the HI ionization rate measurements in Section |3.6| . 
In Section ^ we present Lyman a forest simulations and 
investigate curve-of-growth effects in the proximity effect 
analysis. Section |^ is a summary and discussion of our re- 
sults and the possible systematic effects entering into our 
analysis. 

Throughout this paper, we assume values of 75 km s^^ 
Mpc^^ for Ho and 0.5 for qq. 

2. DATA 
2.1. Spectrophotometry 

Spectrophotometry of 12 sample objects in the spectral 
region between Lya and C IV emission was obtained at the 
Steward Observatory (SO) Bok Telescope with the Boiler 
and Chivens (B&C) Spectrograph and the 1200 x 800 
CCD on the nights of September 22, 1992, November 29, 
1994, and March 28, 1995. Observations were made with a 
400 1 mni^^ grating with Ab=4889 A in the first order and 
a 4.5" slit. Spectrophotometry of the object 1422-f 231 was 
obtained at the SO B&C using a 600 1 mm~^ grating with 
A;,=6681 A in the first order and a 1.5" slit on April 22, 
1996; and the object 1603-1-383 was observed by A.D. as 
part of the Hamburg/CfA Bright Quasar Survey on July 
4, 1995 with the Fred Lawrence Whipple Observatory 1.5- 
meter Tillinghast telescope and FAST spectrograph, using 
a 300 1 mm~^ grating with Af,=4750 in the first order and 
a 3" slit. See Table |l| for a summary. 

All observations except those of 1422-1-231 and 
1603-1-383 were made with the slit set at the parallactic 
angle. This should not seriously effect the spectrophotom- 
etry of 1603-1-383 as it was observed at a small airmass. 
Additionally, however, the observation of 1422-1-231 was 
made with a slit width that is somewhat small for optimal 
spectrophotometry. In any case, as discussed further be- 
low, both 1422-^231 and 16034-383 are excluded from the 
proximity effect analysis due to the fact that 14224-231 is 
a gravitational lens and the presence of associated absorp- 
tion in the spectrum of 16034-383. Any small errors in the 
spectrophotometry of the 74 objects used in the proximity 
effect analysis should not significantly bias the results of 
this work. 

Object spectra were bias corrected and extracted using 
standard IRAF packages using He-Ne-Ar and quartz cali- 
bration exposures taken at each telescope position to per- 
form the wavelength calibration and to correct for pixel- 
to-pixel variations, respectively. The data were then flux 
calibrated using standard star exposures. The column den- 
sity of Galactic neutral hydrogen along the line of sight to 
each object was found using the program GOLDEN, made 
available by J. M'^Dowell; and the spectra were corrected 
for the Galactic reddening calculated from the relation 
Nhi/E{B -V)= 4.8 X 10^1 atoms cm'^ magnitude^! 
(Bohlin et al. 1978). The spectra and the power law con- 
tinuum fits are shown in Figure ^ 

2.2. QSO Systemic Redshifts 

For the present absorption line sample, the QSO narrow 
emission lines discussed above all lie redward of ~7600 A, 
and into the near infrared. Spectra of four objects in this 
sample were obtained at the MMT with the infrared spec- 
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trometer FSpec (Williams et al. 1993) on May 20, 1994 
(1207+399 and 1422+231) and April 1, 1996 (1408+009, 
and 1435+638) using a 75 1 mni~^ grating and a 1.2" 
slit giving a resolution of ^^34 A in the K band. A se- 
ries of exposures of each object was taken. Between each 
exposure, the object was moved along the slit. The to- 
tal integration time is listed in Table |[ One object, 
0836+710, was observed on March 28, 1995 with the B&C, 
the 1200x800 CCD, a 300 1 mm^i grating with Ah==6693 A 
in the first order, and a 4.5" slit. Infrared spectra of eight 
objects in this sample, 0000-263, 0014+813, 0636+680, 
0956+122, 1159+124, 1208+101, 2126-158, were obtained 
using FSpec, OSIRIS on the Cerro Tololo Inter- American 
Observatory 4 m telescope, and CRSP on the Kitt Peak 
National Observatory 4 m telescope (Kuhn 1996). A sum- 
mary of these observations is given in Table g and the 
spectra are displayed in Figure 0. 



3. LYa FOREST STATISTICS FOR Zabs • 
PROXIMITY EFFECT 



Zem- the 



3.1. Spectrophotometry 



In order to perform the proximity effect analysis, the 
flux of each QSO at the Lyman limit is needed. The spec- 
trophotometry data discussed above was used for this pur- 
pose. A power law of the form ^ v~°' was fit to the 
continua of these objects. The straight hue fit to \og{fy) 
vs. \og{v) was done using a robust estimation technique; 
and emission lines found by visually inspecting the spec- 
trum were excluded from the points used in the fit. The 
measured flux at 1450 A and the value of a derived from 
this fit were used to determine the flux at 912 A. For the 
objects we did not observe, we proceed as follows. If a flux 
measurement at a rest UV wavelength other than 912 A 
exists along with a published spectral index, we use these 
to extrapolate to the Lyman limit. If no spectral index 
is available, we use the value of 0.46 (Francis 1996). The 
object 2134+004 has a variable continuum (Perez et al. 
1989, Corbin 1992). Therefore, although we have spec- 
trophotometry from our own observations of this object, 
we take the flux measurement of these authors from their 
averaged spectrum produced from observations made over 
several months. We use this with the spectral index we 
derive to extrapolate to 912 A. 

If no rest UV spectrophotometry of an object exists, we 
estimate at 5500 A (observed) from the V magnitude 
given in Table 1 of Paper I with an extinction correction 
applied. The extinction correction was calculated using 
the column density of neutral hydrogen from GOLDEN 
and the Seaton (1979) re-normalization of the composite 
UV-optical reddening curve of Nandy et al. (1975, and 
references therein). A rest- frame composite QSO spec- 
trum (Zheng et al. 1997) with an arbitrary flux scale 
was redshifted by the appropriate amount for each object. 
The flux in the V fllter was calculated by convolving this 
spectrum with the V filter transmission as a function of 
wavelength. A scaling factor was calculated so that when 
the redshifted QSO composite spectrum was multiplied by 
this factor, the resulting magnitude matched the magni- 
tude listed in Table 1 of Paper I. The flux at 1450 A was 
then taken from this scaled spectrum and this flux was 
extrapolated to the Lyman limit using the spectral index 
given in Table 0. A zero point flux density for the V filter 



of 3.81 X 10^2" ergs s^^ cm"^ jj^-^ (Johnson 1966) was 
used. 

The asterisks in Table || mark QSOs which are known 
lenses or which show associated absorption in their spec- 
tra. Associated absorption is defined to be any Lyman 
a absorption within ~5000 km s~^ of the QSO redshift 
which also shows metal lines. (See Paper I for a descrip- 
tion of the metal line systems identified in each QSO spec- 
trum.) These objects were excluded from the proximity 
effect analysis on the grounds that gas associated with 
the QSO or QSO host galaxy is not part of the general 
intergalactic medium and bulk motions within this gas 
may skew the results. The spectrophotometric properties 
adopted for the 59 QSOs from the literature are listed in 
Table 5 of B94. 

3.2. Number of Lines with Zabs ~ ^em 

The first method we use to demonstrate the proximity 
effect is to compare the number of lines predicted if there 
was no effect from the equation 



— ^Ao{l + zy. 
dz 



(1) 



with the number of lines counted in the spectrum as a 
function of distance from the QSO, 

AJV ^Mpred-J^obs- (2) 

The number of lines predicted is found by integrating 
Equ. |l|. 



pred 



A) 

7+1 



{{l + Z^a.r+' -{l+Zra^nr+'). (3) 



The bins in luminosity distance from the QSO are defined 
according to the relation, 

= 16^^-5 (1+^1)5/2 ^"' Mpc. (4) 

We use h = 0.75. Figure ^ plots the distribution in z and 
Lyman limit luminosity of the QSOs in our sample. 

The dataset was divided into low luminosity and high 
luminosity subsamples at log[L(i'o)]=31.1, such that there 
were equal numbers of objects in each subsample. The 
Lyman limit luminosity of each object was calculated ac- 
cording to the expression 



(5) 

( 1 -|- Zem ) 

where the luminosity distance to the QSO, (11 is given by 
c{qoZ + {q„-l)[il + 2qozy/^-l]} 



(6) 



for qo > 0. In this paper, we use a value of 0.5 for 
go- Figure ^ plots the fractional deficit of lines, {Mpred — 
J^obs) / ■N'predi for the total sample and the high and low 
luminosity subsamples. 

For the total sample, a 5.5cr deficit of lines is found in the 
0-1.5 Mpc bin. The low luminosity subsample shows 
a deficit of lower significance (3.6(t) than the high lumi- 
nosity subsample (4.6(t). These deficits are expected for a 
proximity effect caused by enhanced ionization of HI from 
the QSO flux; and the marginally higher significance for 
high luminosity objects further suggests that this picture 
is legitimate. 
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3.3. Photoionization Model 

We follow the formalism outlined in BDO to calculate 
a value of the mean intensity of the ionizing background 
in the redshift range 1.7 < z < 3.4. The column density 
of a Lya absorber in the immediate vicinity of a QSO will 
be modified from the value that it would have if the QSO 
were not present. The amount by which the column den- 
sity of HI will be reduced due to ionization by UV photons 
from the QSO is given by 



N = No{l+ujy 



(7) 



where A'" is the observed column density of the absorber, 
and Nq is the column density that the absorber would have 
if the QSO were absent. The column density distribution 
of the general Lya absorber population was been shown 
to follow a power law over several orders of magnitude in 
column density, 

MocN-f^, (8) 

which, for a fixed limiting column density, Nthr, (corre- 
sponding to the limiting rest equivalent width) can be 
integrated to give the total number of lines with col- 
umn densities equal to or larger than the limiting value, 

Af{N > Nthr) = Thus, a proximity effect- 

corrected redshift distribution for a fixed rest equivalent 
width threshold can be derived: 



'^^A^[i + zr[i+^(z)] 

az 



(9) 



where w represents a flux-scaled distance of each cloud 
from the QSO 

uj= ^ °; . 10 

47rJ(i/o) 

Here, F'^{i'o) is the Lyman limit flux density due to the 
QSO at the position of a given absorber. 



(11) 



where tl is the luminosity distance between the QSO and 
the absorber. We remove the dominant dependence of the 
line density on redshift by introducing a coevolving coor- 
dinate, Xj, given by 



(1 



Vdz. 



(12) 



If no proximity effect existed, the number of lines per co- 
evolving coordinate would be expressed as 



dJV/dX^ = Aa. 



(13) 



In this analysis, we use a value for /3 of 1.46 from of Hu 
et al. (1995) based upon high S/N, high resolution spec- 
tra of four QSOs at z w 3, consistent with the value of 1.4 
found by Dobrzycki & Bechtold (1996), hereafter DB96, 
from simulations of Lya forest spectra in QSOs at z w 3. 
The value of this parameter is an important factor in the 
ionization model. B94 found that changing the adopted 
value of /3 from 1.7 to 1.4 caused the derived value of J_2i 
to decrease by a factor of ~3. Giallongo et al. (1996) find 



that a double power law provides a better fit to the ob- 
served column density distribution in their high resolution 
spectra than a single power law. The form of their double 
power law consists of a break at Nhi = 10^^ cm^^ and 
values of /3 above and below this break of 1.8 and 1.4 re- 
spectively. For this analysis, however, we will use a single 
power law, as the data of Hu et al. (1995) do not require 
the double power law form. 

The procedure consists of assuming a form for J{vq) as a 
function of z, dividing the lines into the appropriate uj bins, 
and finding the parameters of the assumed form of J(t'o) 
that gives the lowest between the binned data and the 
ionization model. Since no work to date has shown that 
J{vo) evolves significantly with redshift over the range of 
our sample objects, we will treat the case that J(i'o) is 
constant over the redshift range of the data. 

Figure H plots with respect to the constant J(i'o) pho- 
toionization model versus log[J(t'o)] and Figure ^ plots the 
coevolving number density versus u) for the lowest value 
of J{vq) for each subsarnple. The results of this analysis 
are summarized in Table and are discussed in more detail 
in Section ^. 

3.4. Maximum Likelihood Analysis 

In addition to the standard BDO analysis, we also used 
a maximum likelihood method outlined by KF93 to mea- 
sure the extragalactic ionizing background in a manner 
that avoids binning of the data. One constructs a likeli- 
hood function of the form 



L = Y{f{Na,Za)\{exp[ 

a Q 



N'- 



dz I f{N,z)dN], 
(14) 

where the subscripts a and Q refer to absorbers and QSOs 
and where f{N,z) is the standard equation for the dis- 
tribution of Lyman a absorbers in column density and 
redshift. 



f{N, z) = AN-f^il + zy[l + Lo{z)]-^'^-^\ 



(15) 



The parameter lo is defined as above, but here, the 
normalization in terms of in Equation |l| is given by 
Ao{Num/NQ)P~^{l3-iy^ . With the exception of the case 
in which a variable threshold is used, Nmin for each QSO 
is the column density which, according to the curve-of- 
growth adopted (see KF93), corresponds to an equivalent 
width of 0.32 A, 2.62 X lO^^^ cm^^. 

Instead of using the method outlined by KF93 whereby 
the parameters A,/3,7, and J{vq) are all found by mini- 
mizing -ln(L) where L is given by the likelihood function 
above, we chose to take the parameter 7 from a separate 
maximum likelihood solution to Equation]^ (sec Paper I.) 
Since our spectra are more highly blended than the low 
redshift data used by KF93, we choose not to determine 
(3 directly from our data using line equivalent widths and 
the curve-of-growth and instead adopt a value found from 
high resolution spectra. As described in the previous sec- 
tion, we take (3 to be 1.46 (Hu et al. (1995) and solve 
for A by requiring f(N,z) to give the observed number of 
lines in the regions of the QSO spectra unaffected by the 
proximity effect. 

We ran two tests on this set of algorithms. The first of 
these was to attempt to reproduce the results of KF93 with 
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the dataset they used from BahcaU et al. (1993) . Next, 
we used a high redshift subsample of our complete dataset, 
the DB96 sample, to compare the results of the maximum 
likelihood analysis and the EDO analysis to each other and 
to independent checks on these values (B94, Giallongo et 
al. 1996). 

We were able to reproduce the results of KF93. Using 
their Sample 2, the BahcaU et al. (1993) sample minus 
one BAL QSO, PG 0043+039, we obtain (7, (3, log(A)) = 
(0.23, 1.47, 7.74) and log[J(z^o)]=-23.0t|^;6 for 6 = 35 km 
s^^. These agree with the values they find, (7, /?, log(A))= 
(0.21, 1.48, 7.74), and the errors in these values, ~ 0.06 
f7/3 ~ 0.05, and criog{A) ^ 0.1. Their result for log[J(fo)] 
for this sample is -23.3lg'5. 

The high redshift subsample we created consisted of 
518 lines from the 15 objects from DB96 that do not 
show associated absorption. The QSOs have redshifts 
between 2.52 and 3.38. Using our maximum likelihood 
program to solve for the Lya forest statistics, we find 
7 = 1.926 ± 0.656, and \og{A) = 7.03 for 7V„„,=2.6 x 
10^"* cm~^ and (3= 1.46. This subsample does give similar 
results in the BDO and the maximum likelihood cases, 
log[J{i^o)]=-2lAQtlfg and log[j(z.o)]=-21.58l°jo, re- 
spectively. (See rows 1 and 2 of Table 0.) These val- 
ues agree well with the Giallongo et al. (1996) result of 
log[J(i/o)]=-21.30 ± 0.7 for z=l. 7-4.1. 

The software we used for the maximum likelihood anal- 
ysis uses all regions of the QSO spectra between Zmin, 
specified by the spectral coverage or by Ly/3 emission, and 
Zmax , specified by Lya emission. Though it does not count 
lines associated with identified metal line systems, it does 
not exclude the regions of the spectrum where these lines 
lie. To ensure that this does not have a significant effect 
on our resultant solution for the background, we tested a 
program that does exclude regions of the spectra in the 
same way that our BDO-style software does. The change 
in the result was indeed insignificant; but taking these ex- 
cluded spectral regions into account and binning the data 
in the same way the BDO-style software does brings the 
maximum likelihood and the BDO method results into ex- 
cellent agreement. 

Figure |7| plots the log of the ratio of the likelihood func- 
tion to the maximum value versus log[J(j/o)]; and Figure |^ 
plots the coevolving number distribution of Lya lines with 
respect to u just as in Figure ^. The results of this anal- 
ysis are also summarized in Table ^ and discussed further 
in Section 

3.5. Systemic QSO Redshifts 

One of the major uncertainties in the proximity effect 
analysis is in the systemic redshifts of the QSOs. If the 
true redshift of a QSO is higher than the value used in the 
analysis, any given cloud is further away from the QSO 
than assumed. Hence, the influence of the QSO at this 
cloud is less than inferred and the value of J{i'o) in reality 
is lower than the one derived. 

For the data presented in Figure ^, an average of sev- 
eral cursor settings at the peak of the emission line was 
used to determine the line centers. More detailed fits were 
not done as our purpose lies mainly in determining if any 
gross shifts between Lya and the Balmer lines/ [OIII]/Mg 
II exist for our data; but we found no significant difference 



between this method and making Gaussian fits to the up- 
per 50% of the emission line profiles. 

Lya redshifts were measured from the absorption line 
spectra when the entire Lya profile was observed, in the 
same way as was done for the Balmer, [OIII], and Mg II 
lines. Table || lists the adopted best redshift value for each 
emission line for each object supplementing our measure- 
ments with measurements from the literature. 

Laor et al. (1994) and Laor et al. (1995) found, from 
a sample of 13 QSO spectra from the Faint Object Spec- 
trograph on Hubble Space Telescope between redshifts of 
z ~ 0.16 and z ~ 2.0, average velocity shifts between 
[OIII] A5007 and Lya, Mg H, and H/3 of 200± 150 km s-\ 
—85 ± 130 km s~^, and —75 ± 110 km s~^, respectively. 
This agrees with the Corbin & Boroson (1996) result for 
48 objects with 0.03 < z < 0.77. They found mean [OIII]- 
Lya and [0III]-H/5 shifts of 191 ±101 km s~i and -75 ±57 
km s^^. Thus, Lya is blueshiftcd with respect to [OIII] 
by ~200 km s^^, while Mg II and H/3 are marginally red- 
shifted with respect to [OIII]. Tytler & Fan (1992) find a 
mean [0III]-H/3 shift of -15 ± 37 km s'^ from 8 QSOs 
with redshifts between ~0.3 and ~0.6 and conclude that 
both Balmer lines and narrow forbidden lines give red- 
shifts within 100 km s~^ or less of the QSO systemic red- 
shift. They then find a blueshift of Mg II with respect to 
[0III]/H/3 for 100 QSOs of 101 ± 47 km s'^ which they 
use as a secondary systemic redshift zero point in their 
analysis of a large QSO sample. The magnitude of the 
blueshift of Lya with respect to [OIII] /H/3 that they de- 
rive IS 172 ± 17 km s~^ The data of Nishihara et al. 
(1997) for five QSOs at z ^ 1.5 show a negligible redshift 
of Mg II with respect to [OIII], 31 ± 411 km s"^ However 
these five objects show a somewhat larger redshift of H/? 
with respect to [OIII] A5007, equalling 260 ± 522 km s-\ 
consistent with the fact that these objects have high lu- 
minosities. M'^Intosh et al. (1999b) use the near-infrared 
spectra of QSOs at 2.0 ^ z ^ 2.5 presented in M'^Intosh 
et al. (1999a) to examine the redshift differences between 
[OIII] and H/3. They supplement their data with data from 
the literature to measure the redshift differences between 
[OIII] and Mg II. They find that on average, H/3 is red- 
shifted relative to [OIII] by 520 ± 80 km s'^ for 21 of their 
sample objects, while Mg II lies within 50 km s^^ of the 
redshift of [OIII] for 12 sample objects. 

For our sample, we find that Lya is blueshiftcd with re- 
spect to [OIII] A5007 by 382 ± 1160 km s'^ for 19 QSOs. 
Mg II emission is blueshiftcd by an average of 338 ±901 
km with respect to [OIII] on the basis of seven mea- 
surements. We find that H/3 is redshifted by 642 ± 740 
km s^^ with respect to [OIII] on the basis of five measure- 
ments. Including three Ha redshifts listed in Table ^ with 
these H/3 redshifts leads to a 507 ± 615 km s~^ redshift 
of Balmer lines with respect to [OIII]. This shift is larger 
than that discussed above for low reshift QSOs. However, 
it is consistent with the Nishihara et al. (1997) H/3 shift 
for high luminosity QSOs. Combining our data with that 
of these authors, we find that Mg II is blueshiftcd with re- 
spect to [OIII] by 184 ± 735 km s^^. Including the data of 
M'^Intosh et al. (1999b) that is not already in our sample 
gives a blueshift of 95 ± 603 km s~^. Similarly, combin- 
ing our data with that of Nishihara et al. (1997), we find 
that H/3 is redshifted with respect to [OIII] by 451 ± 636 
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km s^^. After supplementing this combined data set with 
the data of M'^Intosh et al. (1999b), the redshift becomes 
379±516 km s~^. Finally, combining the data of M'^Intosh 
et al. (1999b) with ours gives a Lya blucshift of 418 ± 920 
km s~^ with respect to [OIII]. 

As has been noted in previous work, the standard er- 
ror in the mean velocity shifts is quite large, on the order 
of or exceeding the value of the shift itself. We estimate 
that the wavelength calibration errors in our data con- 
tribute a ~10-30 km error in the derived redshifts; 
and the spread in different redshift measurements of the 
same species (e.g. Balmer lines or [OIII] A4959 and A5007) 
for the same object is typically 100-200 km s~^. The ob- 
served spreads in the velocity differences of the Lya, Mg 
II, and Balmer emission lines with respect to the QSO sys- 
temic redshifts are much larger than this, indicating that 
it is intrinsic to the QSO population. Figure ^ shows his- 
tograms of the emission line redshift differences between 
[OIII] and Lya, [OIII] and Mg II, and [OIII] and Balmer 
lines. Our results are plotted with those of Laor et al. 
(1995) and of Nishihara et al. (1997). Our sample shows 
no well-defined mean [OIII]-Balmer shift, just a large scat- 
ter in the measurements included. Our sample also shows 
a large range of [OIII]-Lya and [OIII]-Mg II shifts with 
no well-defined mean value. Nonetheless, the mean trend 
is that the [OIII]-Lya shift is different from zero by lAa 
for our data, less than the 3.5cr significance found by Laor 
et al. (1995). The [OIII]-Balmer line shifts for both our 
data set and for our data combined with that of Nishihara 
et al. (1997) are more significant, 2.7cr and 2.8(7 respec- 
tively. The [OIII]-Mg II shift is consistent with zero in a 
mean sense, but with large scatter. Thus, though better 
statistics are desirable, it seems that for these high red- 
shift and relatively high luminosity objects, Balmer lines 
are not good indicators of the QSO systemic redshift. For 
the purposes of this study therefore, we treat only the red- 
shifts found from [OIII] A5007 for 19 objects in our sample 
and Mg II for 16 objects in our sample as systemic QSO 
redshifts. 

3.6. The HI Ionization Rate 

The HI ionization rate due to a source of UV flux is 
formally given by the equation: 



r = 



hv 



dv s ^. 



(16) 



The calculations of the mean intensity of the ionizing back- 
ground to date have made a critical assumption, namely 
that the spectrum of the background and the spectra of 
the individual QSOs are identical. This allows the ex- 
pression w = 

YQ/T^9 to reduce to the ratio of the Ly- 
man limit flux density of the QSO, J^{ua), to that of 
the background, J^^{vq), for each line (BDO). Since the 
IGM reprocesses the radiation emitted from QSOs, this is 
not strictly true (Miralda-Escude & Ostriker 1990, Madau 
1991,1992, Meiksin & Madau 1993, Haardt & Madau 1996, 
Fardal et al. 1998). Furthermore, the value of P''^ is of 
particular interest as it can be used to infer the value 
of rif, by comparing the distribution of flux decrements 
in high resolution QSO spectra to Lyman a forest sim- 
ulations (Ranch et al. 1997). Therefore, we repeat the 
standard BDO analysis without making this assumption. 



ie. using u — T'^ jT^s and solving for the HI ionization 
rate from the metagalactic background radiation. The 
ionization rate for each QSO was calculated using Equa- 
tion H, where aHi{v) = 6.3 x cm^ and where 
jQ{v) = jQ{j,f)){±-y°'. For each QSO, J'^ii'o) is the 
same value used in the standard analysis used to solve for 
j''^{i^o), and a is given in Table ||. For some objects, no 
a listed in this ta ble a nd a value of 0.46 was used, as de- 
scribed in Section 3.1. As before, the best value will be 



the one that gives the lowest between the model with 
/3=1.46 and the binned data. We use the narrow line red- 
shifts for each QSO discussed above and add 400 km s""'^ to 
each QSO redshift measured from the Lyman a emission 
line. 

Haardt & Madau (1996) present a Gaussian fit to their 
model for the evolution of F with redshift, 



A{l + zfexp[~{z~Zcf/S] 



(17) 



that agrees with their detailed model for the background 
to within 10% over the range < z < 5. The best fit 
parameters they derive for 50=0.5 are A— 6.7 x 10"^^ s~^, 
B=0A3, Zc=2.30, and 5=1.95. Fardal et al. (1998) fit 
their model for the background with the parameter sets 
^=5.6 X 10-" s-i, B=0.60, Zc=2.22 , and 5=1.90 and 
^=1.26 X 10-12 s-\ 5=0.58, Zc=2.77 , and 5=2.38 for 
the Ql and Q2 luminosity functions, of Pei (1995) respec- 
tively. Incorporating this expression for T(z) with these 
three different sets of parameters into the BDO style anal- 
ysis allows us to determine which of these models fits our 
data best. The results are listed in Table |6| and are dis- 
cussed in greater depth below in Section 

4. SIMULATIONS AND THE CURVE OF GROWTH 

Simulated Lyman a forest spectra for the DB96 sample 
only were produced using the software described in that 
paper. The simulation input 7 was changed slightly to re- 
flect the maximum likelihood value found by the software 
used in the analysis described in Paper I. The normaliza- 
tion was chosen to give matching amounts of total absorp- 
tion in the real and simulated spectra. The parameters 
used were 7 = 2.069, Aq = 4.835, /? = 1.46, log(NHi„,.J = 
13.0, log(NHi„„J = 16.0, <b>= 28.0 km s^^, ai,= 10.0 
km 3-1, and 6cut= 20.0 km s^^. 

The proximity effect was included in these simulations 
by simply modifying each cloud's column density accord- 
ing to equations |^ and The value of log[J(i^o)] from 
the BDO type analysis on the DB96 sample is -21.40^|!) gg. 
Values of -19.0, -20.0, -21.3, -22.0, and -23.0 for log[J(^'o)] 
were input and the analyses described above were used to 
recover that J{vo). Two examples of the simulated spectra 
are shown in Figure |l^. 

The analysis considers all lines above a fixed equivalent 
width threshold of 0.32 A. Thus, as the column densities 
of lines are modified by the QSO flux from their expected 
values in the absence of the proximity effect, the equivalent 
widths of the lines will change according to the curve-of- 
growth. If a line is saturated, changing its column density 
will have little effect on its equivalent width, since it lies on 
the flat part of the curve-of-growth where W oc y^log{N). 
This will mean that for a given equivalent width cutoff in 
the data, this line will not drop out of the sample as the 
proximity effect is turned on in the simulations. Since the 
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line deficit will be less than expected for a given input value 
of J(fo), the proximity effect will appear less pronounced 
and the true J(z/o) will be overestimated. We found this 
to be the case from our simulations. As Figure |l^ illus- 
trates and Table ^ summarizes, though the values of J{vo) 
recovered from the simulated data were usually consistent 
with the input values within the Icr confidence limits, they 
were systematically larger than the input values by up to 
a factor of 3. The largest input values of log[J(i'o)], -19.0 
and -20.0, give the largest discrepancy between this input 
value and the log[J(j/o)] recovered from the BDO analysis 
performed on the simulated spectra. The smallest input 
value of log[J(i'o)], -23.0, gives the smallest discrepancy 
between the input and recovered values. However, the 
la confidence limits on this fit are also relatively small, 
making it the only trial which does not recover the input 
log[J(i/o)] to within those limits. 

To demonstrate the effect. Figure |l^ compares the simu- 
lated line equivalent widths with and without the proxim- 
ity effect included. The column density of each line from 
the simulated spectra line lists with no proximity effect 
were modified according to equations I and |lO|. Fi gures 
12(a-e) plot the nonproximity effect rest equivalent width 
Wno-PE versus the ratio of the proximity effect and non- 
proximity effect equivalent widths, WpE /Wno-PE- The 
solid line delineates the detection threshold for the lines in 
the list for which the proximity effect is included, WpE^ 
0.32 A. Absorption lines that fall above this line were not 
removed from the sample when the proximity effect was 
turned on, while those below it disappeared. For a given 
set of QSOs with fixed Lyman limit lumosities, such as this 
one, the proximity effect signature in their spectra will 
become less pronounced as the ambient UV background 
increases. Therefore, as log[J(i'o)] increases from -23.0 to 
-19.0, the magnitude of the proximity effect decreases, and 
the proximity effect line list differs less and less from the 
nonproximity effect line list. 

5. RESULTS AND DISCUSSION 

Table ^ lists the best fit values of J(i'o) found for var- 
ious subsamples of this dataset using both the canonical 
BDO and the maximum likelihood methods. For the BDO 
method, the Icr confidence limits are found from a Ax^ 
of 8.18 for 7 degrees of freedom. The maximum likeli- 
hood method la confidence limits derive from the fact 
that \n{L/Lmax) is distributed as The total sam- 

ple consisting of 74 QSOs with all QSO redshifts based on 
the Lya emission line gives a best fit value of log[J(t'o)] of 
-20.90155^8 for the BDO analysis and -20.83l5J;2o for the 
maximum likelihood analysis. 

As the results in Table demonstrate, using narrow line 
redshifts for 35 of the 74 QSOs for which they have been 
directly measured and Lya redshifts for the rest does not 
change the result. However, when 400 km s^^ is added to 
the Lya redshifts of the objects with no measured_narrow 
line redshift, a value for log[J(i'o)] of 
rived using the BDO method and \og[J{h'o 
is found using the maximum likelihood method. Recall 
that the mean blueshift of Lya with respect to [OIII] for 
the 19 objects in this paper with [OIII] A5007 measure- 
ments was found to be '^400 km s~^. This decrease in 
the mean intensity of the background derived when larger 



-21.15tj^ 4^ is de- 
-21 1 7+0- 19 



QSO redshifts are used is to be expected, (cf. Section 3.5 ) 
Because this measurement of the background accounts for 
the systematic blueshift of the Lya emission line with re- 
spect to the systemic redshift of each QSO, we consider 
it to be our best estimate for the mean intensity of the 
background at the Lyman limit. 

These measurements have been made, however, using a 
photoionization model with somewhat unrealistic assump- 
tions, particularly that Lya absorbers are isothermal and 
are composed of pure hydrogen. For clouds with a primor- 
dial He abundance and which are in thermal and ionization 
equilibrium, Using CLOUDY to model the ionization state 
of absorbers with a metal adundance of 10^^ solar (Cowie 
et al. 1995, Tytler & Fan 1994) as a function of a;, we find 
that the neutral fraction, x, is proportional to (l + uj)^^''^^ . 
This implies that 

^ = Aoii + zr[i+io{z)r'-''(^-'^ (18) 

In this scenario, the optimal value found for log[J(i/o)] is 
-21.101q28. This value is marginally larger than the value 
discussed above, found under the assumption of absorbers 
composed of pure hydrogen; but it is not significantly dif- 
ferent, so we conclude that the absence of metals in the 
BDO model has not drastically affected our measurement 
of the background. 

It is worth noting that 16 objects in our sample of 
objects with no associated absorption show evidence for 
damped Lya absorption: 0058-1-019, 0100-1-130, 0334-204, 
0913-1-072, 0938-hll9, 0952-^338, 0955-^472, 1009-^299, 
1017+280, 1215-1-333, 1247-f267, 1548-1-092, 1946-1-770, 
2126-158, 2233-1-131, and 2320+079. The dust in these 
systems could cause the intrinsic QSO fluxes to be under- 
estimated. This in turn can cause log[J(i'o)] to be un- 
derestimated by up to a factor of 3, in addition to the 
sources of error discussed above (Srianand & Khare 1996). 
Only six of these objects, 0334-204, 0938+119, 0955+472, 
1215+333, 2126-158, and 2233+131, appear in our low 
luminosity subsample, suggesting that this subsample is 
not preferentially heavily dust-obscured. Nevertheless, the 
BDO analysis was performed on all 16 objects exhibiting 
damped Lya systems; and found the best fit value for 
log[J(i'o)] to be -21A5toi°3, ^ factor of 1.9^^;^ lower than 
the value obtained for the sample as a whole. This does 
not allow us to say anything significant about the presence 
or absence of dust, so we will neglect its influence. 

Dividing our line sample into subsamples of high (z > 
2.5) and low (z < 2.5) redshift lines, we find marginal 
evidence for evolution in the intensity of the background, 
namely that the maximum likelihood background inten- 
sity is lower by a factor of about 1.91^^ 4 lower redshift. 
The BDO results corroborate this, but with larger uncer- 
tainties. The factor by which J{i'o) is found to be lower 
at lower redshifts is 2.5l2^2^- Gravitational lensing could 
mimic a trend with redshift with about the same order 
of magnitude, if the high redshift subsample contains a 
significant number of unknown lenses. However, Figure |^ 
suggests little if any trend for high luminosity objects to 
exist at high redshifts in our sample; and the results of 
Section 3.2 indicate that the high luminosity objects do 
show a somewhat stronger proximity effect despite the fact 
that the measured background at high redshift appears to 
be higher. No other studies have found this evidence of 
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redshift evolution in the background, so we regard it as 
tentative; and note that it will be interesting to see in 
future work if this trend can be shown to be real and if 
it extends smoothly to the low values of J{vo) found at 
redshifts less than 1.5. 

Since we find high luminosity objects do not exist pref- 
erentially at high redshift in our sample, a simple test can 
be done to determine whether or not there is a signifi- 
cant number of lensed objects in our sample. If the high 
luminosity QSOs are indeed intrinsically more luminous, 
and the proximity effect is a purely photoionization-driven 
phenomenon, these objects should show a m ore prominent 
proximity effect. The results of Section 3.2 suggest this is 
the case. However, in the analysis, this larger line deficit is 
normalized to the higher Lyman limit luminosities of this 
subsample. Therefore, one expects these objects, when an- 
alyzed as a separate subsample, to yield a value of J(i'o) 
that is consistent with that found for low luminosity ob- 
jects if the values of the QSO fluxes are not in error due to 
lensing. If the high luminosity QSOs, or a subset of them, 
are lensed objects, then they are not necessarily intrinsi- 
cally more luminous than the low luminosity QSOs. In this 
case, the influence of the lensed objects on the surround- 
ing IGM will be overestimated and given the observed line 
deficit, the background will also be overestimated. Table ^ 
lists the results obtained for the high and low luminos- 
ity subsamples of our data set. The values obtained for 
these subsamples are equal within the uncertainties. This 
is consistent with there being no significant effects from 
gravitational lensing in our sample. 

5.1. HI Ionization Rate 

We tested a range of values for F, the HI ionization 
rate, using our data. The constant value found to fit the 
data the best is l.9t\:l x 10"^^ This value IS m 

reasonable agreement with that predicted by the QSO- 
dominated model of Haardt & Madau (1996) at this red- 
shift, 1.0 X 10-12 s-i (go=0.5). Using Equation |l| and 
jQ(jy) — and assuming global QSO spec- 

tral indicies of 0, 1.5, and 2, the ionization rate found 
from our data corresponds to log[J(t^o)]= -21.34, -21.17, 
and -21.12, respectively. 

The parameter set (A, B, Zc, (S) found to give the best 
fit to the data is that of Fardal et al. (1998) for the Q2 
luminosity function (1.2 x lO'^^ g-i^ o.58, 2.38, 2.77) 
which, for a redshift of 2.9 yields an ionization rate of 2.7 
X IQ-i^ s^i, in good agreement with our solution, and 
within a factor of ^^3 of the Haardt & Madau result. Thus, 
we conclude that a significant contribution to the ionizing 
background from stellar UV emission is not required at 
this redshift. 

5.2. Curve- of- Growth and Other Systematics 

On the basis of a curve-of-growth argument, one might 
expect that weak lines would show a more prominent prox- 
imity effect than strong lines. We have compared the re- 
sults obtained for a constant equivalent width threshold of 
0.32 A with that obtained for lines with 0.16 A < W < 
0.32 A. Instead of finding a more pronounced proximity 
effect for the weak lines, we find a less significant deficit of 
lines within 1.5 h~^ Mpc of the QSOs. This deficit is 4. Oct, 
versus 5.5ct for hues with W > 0.32 A. As Table | Usts, 



the value of log[J(i^o)] recovered from these weak lines is 
correspondingly higher than that found using strong lines, 
-20A5tofQ versus -21.15t[!i^. Cooke et al. (1997) point 
out that this could be the result of a higher degree of blend- 
ing of weaker lines compared to strong ones in crowded 
spectral regions. The background flux measurement will 
be an overestimate because blending will cause fewer indi- 
vidual lines to be resolved further from the QSO. Because 
the reduction in line density near the QSO will work to 
reduce line blending, the overall effect of line blending will 
be to suppress the true magnitude of the proximity effect 
causing J{vo) to be overestimated, by a factor of 4.5 in 
this case. It is difficult to ascertain whether this effect is 
as strong for lines with > 0.32 A or whether the curve- 
of-growth effect discussed in Section ^ which also causes 
J{vo) to be overestimated, is more important. We expect 
that for lines with W > 0.32 A, the effects of blending are 
reduced somewhat, while the curve-of-growth effects will 
remain a factor. 

We have addressed many of the systematics which could 
possibly have affected our analysis. A treatment of the 
QSO systemic redshifts was integrated directly into our 
analysis and was found to influence the J(i'o) found by 
up to a factor of ~2. Other effects, such as the influences 
of metals and dust, which can cause Ji^o) to be underes- 
timated, and the influences of lensing, line blending, and 
curve-of-growth effects, which can cause J(i^o) to be over- 
estimated, were treated after the fact in an attempt to 
understand the magnitude of their effects on the value of 
J(t'o) derived. The CLOUDY simulations discussed above 
indicate that allowing for an absorber metal abundance 
of 10^2 solar has little effect on the value of J(i'o) found 
from the data. Dust in intervening absorption systems 
may have affected our result. Though we were unable to 
quantify this effect with high confidence, it could be on 
the order of a factor of 2. We assert that QSO flux am- 
plification due to lensing has not significantly biassed our 
result; and we attempt to minimize the effect of blend- 
ing discussed above by using only lines with W > 0.32 A. 
Our result may be susceptible to the curve-of-growth ef- 
fect we addressed through the simulations in Section ^. In 
those simulations, we found that the discrepancy between 
in the input and recovered values of J(i'o) depended up- 
ong the input value of J(i^o) itself. The magnitude of the 
discrepancy corresponding to the J{vq) we found from the 
data was a factor of ^2. We therefore suspect that if our 
result, log[J(i'o)]=-21.15to 43, is systematically biased in 
any way, it is an overestimate of the true background and 
could be in error by up to a factor of 2; though this could 
be balanced somewhat by systematic error due to dust, 
which works in the opposite direction. 

5.3. Comparison with Previous Measurements 

Our value for J(j^o) agrees well with other measurements 
at similar redshift, with the exception of those of B94 and 
Fernandez-Soto et al. (1995) who both derive values four 
times larger than our best value for J(i'o), ~3 x 10-^^ ergs 
cm~2 Hz"-'^ sr-i. The measurement of B94 does not 
take into account QSO systemic redshifts, but she notes 
that if they are blueshifted with respect to Lya by 1000 
km s^i, this would lower the derived value of J(i'o) by 
a factor of 3, bringing it into reasonable agreement with 
our result. The Fernandez-Soto et al. (1995) value is de- 



9 



rived from 3 QSO spectra showing a proximity effect due 
to foreground QSOs. These authors are not able to place 
an upper limit on their measurement, but our value of 7.0 
X 10^^^ ergs s~^ cm^^ Hz^^ sr~^ for J(;^o) is consistent 
with their lower limit of 1.6 x 10~^^ ergs s""'^ cm~^ Hz^^ 
sr~^. In fact, when these authors examine the proximity 
effect in a single QSO spectrum due to the background z^2 
QSO itself, they derive a value for J{vq) of T.of^^g x 10"^^ 
ergs s""'^ cm~^ Hz~^ sr~^, which brings their estimate into 
better agreement with our values for our total sample and 
for our low redshift subsample within their large errors. 
Direct measurements of the background at redshifts ~3- 
3.5 have been made using long-slit spectroscopy of fields 
containing optically thick Lya absorbers in efforts to de- 
tect fluorescent emission the absorbers produce from the 
ionizing radiation field incident upon them (Lowenthal et 
al. 1990, Martinez-Gonzalez et al. 1995). Recent Keck 
telescope observations by Bunker et al. (1998) at 2.5 < z 
< 4.1 have achieved a factor of 2-10 higher sensitivity and 
place a firmer direct limit on the background than pre- 
vious work. Their null signal in a 90-minute integration 
with a 3' slit sets an upper limit on J{i'o) of 2 x 10 
ergs s~^ cm~^ Hz~^ sr~^. 

Cooke et al. (1997) claim that the value for the back 
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ground at z~4 is between their value of 8.0 



+8.0 
-4.0 



X 10- 



-22 



ergs s ^ cm ^ Hz ^ sr ^ and that of WiUiger et al. (1994), 
1.0-3.0 X 10^22 g^gg g-i pjj^-2 jj^-i sr-i. 



Our best value 



Hz 



-1 



sr 



of J(i/o) at Z--3, 7.0 X 10"^^ ergs s^^ ci 
is in agreement with this, although within the uncertainty 
there is an allowance for the background to decrease as z 
approaches 4. 

Table 1^ lists these various measurements of Jiiy'o) in the 
literature as well as the Kulkarni & Fall (1993) measure- 
ment at z~0.5. Figure ^ also summarizes the literature 
measurements of J{vo) from Z'--^0.5 to z=4.5. 

The solid curves in Figure |l3| delineate the evolution of 
the mean background intensity as a function of redshift for 
global background source spectral indicies between and 
2, derived from the Haardt & Madau (1996) model for the 
HI photoionization rate as a function of redshift discussed 
in Section Over 90% of our sample QSO redshifts 

lie within the FWHM of the Gaussian in the Haardt & 
Madau (1996) expression using their best fit parameters. 
At these redshifts, the Haardt & Madau (1996) curves in 
Figure |l3| are turning over. Nonetheless, for comparison 
with previous work (B94 and references therein) , we inves- 
tigate a power law redshift dependence of the background 
intensity: 

J{vQ,z) = J{v^,Q){\ + zy. (19) 

Using the EDO method, we executed a crude grid search 
in an attempt to constrain the power law index and nor- 
malization of this power law. The lowest (3.86) between 
the binned data and the EDO photoionization model for a 
power law background was achieved by (j, log[J(;/n, 0)])= 
(5.12, -23.97), shown by a dashed line in Figure 14. Ex- 
tending this solution to low redshift gives log[J(i'o), 0.5)]=- 
23.0, in good agreement with the measurement of Kulka- 
rni & Fall (1993). The solution (j, log[J(i^o, 0)])= (-4.16, 
-18.76) gives the next lowest (4.91); and though it also 
implies mean background intensities over four orders of 
magnitude too high at low redshift, it traces the Haardt 
& Madau model at high redshift, giving log[J(i'o, 4.5)]=- 



21.8, in agreement with the Willigher et al. (1994) mea- 
surement. It is also shown by a dashed line in Figure |4[ 
Fitting parabolas to the regions near the minima in 
both j and log[J(i/o, 0)] gives the error in each param- 
eter for both of these solutions, (5.12 ± 1.96, -23.97 ± 
1.07) and (-4.16 ± 2.36, -18.76 ± 1.31). E94 found a sim- 
ilarly large range of acceptable solutions: — 7 < j < 4 
and —16.5 < log[J(i'o, 0)] < —23.0. The large error bars 
on these fits indicate that the power law fit to the data 
is not well-constrained, due possibly to the fact that the 
mean intensity of the background is turning over at the 
redshifts of our sample objects, as the Haardt & Madau 
(1996) model predicts. 

5.4. Comparison with Models for the Background 

Recent models of the ionizing background include not 
only the integrated emission from QSOs but also a variety 
of other physical processes such as star formation in young, 
high redshift galaxies and attenuation of UV photons by 
Lyof absorbers and Lyman limit systems (Miralda-Escude 
& Ostriker 1990, Madau 1991, 1992, Meiksin & Madau 
1993, Haardt & Madau 1996,Fardal et al. 1998). Madau 
& ShuU (1996) find that the production of metals in Lya 
absorbers may also be a significant contributor to the UV 
background at z « 3. Their contribution may be up to 5 
X 10~^^ ergs s~^ cm~^ Hz~^ sr~^, assuming that the bulk 
of the metals in the Lyman a forest did not form at z >> 
3, and assuming a Lyman continuum escape fraction, /esc, 
from a galaxy of ^0.25. They note, however, that f^sc is 
essentially unconstrained. 

Past debate about how the space density of QSOs 
evolves at high redshift (Koo & Kron 1988; Eoyle et al. 
1991; Irwin et al. 1991; Schmidt et al. 1991; Warren et 
al. 1994; Kennefick et al. 1995) has been clarified by re- 
cent radio surveys (Hook et al. 1995, 1998; Shaver et al. 
1996). This work has demonstrated that the space density 
of radio-loud QSOs decreases rapidly with redshift beyond 
z^3. Since these surveys are unaffected by any presence 
of dust in the intervening IGM; and since they confirm the 
behavior seen in optically selected surveys, they indicate 
that the QSO population is truly declining at high redshift. 
Nevertheless, the discovery of QSOs with redshifts greater 
than 4 has brought better agreement between the values 
of J{vq) found via the proximity effect and the values pre- 
dicted by the models with QSOs primarily contributing 
to the background (Madau 1992, Meiksin & Madau 1993, 
Haardt & Madau 1996). 

Madau (1992) and Meiksin & Madau (1993) estimate 
the QSO UV background by integrating the QSO luminos- 
ity function (Eoyle 1991) and including the effects of atten- 
uation by hydrogen in the IGM. Their estimates however, 
1-3 X 10~^^ ergs s~^ cm~^ Hz~^ sr~^, are still somewhat 
lower than the values derived in this paper. The analysis 
of Haardt & Madau (1996) takes into account the effects 
of various atomic processes leading to the production of 
hydrogen-ionizing photons within Lya absorbers and Ly- 
man limits systems themselves. They conclude that ob- 
served QSOs can account for number of ionizing photons 
required by the proximity effect at z<^4. These authors 
find a value of log[J(i'o)] equal to ~-21.4 at z=3, in good 
agreement with the value found in this paper at similar 
redshifts. The solid lines in Figure 13|show the results from 
the Haardt & Madau (1996) model for two different values 
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of the global background source spectral index. The lower 
and upper curves show the evolution of the background for 
indicies of and 2 respectively. The literature measure- 
ments at redshifts between 1.7 and 3.6 agree well with the 
model predictions. The Z'-~^0.5 measurement of Kulkarni & 
Fall (1993) falls below both model curves and the z=4.5 
measurement of Williger et al. (1994) falls above them. 

Madau, Haardt, & Rees (1999) revisit the issue of the 
contribution of high redshift, star-forming galaxies to the 
ionizing background in light of recent work identifying such 
objects at 2 < z < 4. (Steidel et al. 1996a,b; Madau et al. 
1996; Lowenthal et al. 1997) They calculate the critical 
photoionization rate necessary to reionize a non-uniform 
intergalactic medium as a function of redshift. This is 
compared to the expected contributions from QSOs and 
young, star-forming galaxies. There are uncertainties in 
estimating both of these. The QSO luminosity function 
at z > 4 must be extrapolated from that at lower red- 
shifts. There is also still some debate between theory and 
observations, eg. of the Hubble Deep Field, on the sub- 
ject of a population of low-luminosity QSOs (see Madau 
et al. 1999 and references therein) which could cause the 
QSO luminosity function to steepen with lookback time, 
making up for the dearth of observed objects at z > 4. 
The estimation of the galaxy contribution of ionizing pho- 
tons is limited by poor knowledge of luminosity function 
of Lyman-break galaxies at z > 4 as well as by the lack 
of constraints upon fpsc- Nevertheless, the results are in- 
triguing. Assuming that /esc=0.5, Madau et al. (1999) 
find that the contribution of hydrogen-ionizing photons 
from star-forming galaxies z ~ 3 could exceed that from 
QSOs by a factor of more than 3. However, the QSO con- 
tribution at this redshift is sufficient, according to these 
estimates, to ionize the IGM at this redshift. Deharveng 
et al. (1997) estimate a much lower fesc at z=0, less than 
1%, based on the local galaxy Ha luminosity density. Fur- 
thermore, Devriendt et al. (1998) make an independent 
estimation of the galaxy contribution to J(t'o) assuming 
damped Lya systems to be the progenitors of present day 
galaxies. Their semi-analytic models include a treatment 
of not only HI absorption of Lyman limit photons in the 
intervening IGM, but also of HI and dust absorption in 
the interstellar medium of the photon-producing galaxies. 
Their results show that constraining /esc in this way yields 
a much lower contribution to the UV background from 
galaxies at z > 2. At z ^ 2.5, their estimated QSO con- 
tribution to J(t'o) is 3 orders of magnitude greater than 
that expected from galaxies. Our measurement of J{vo) is 
consistent with the UV background being QSO dominated 
in the models of both these authors and Haardt & Madau 
(1996). 

In the models of Madau et al. (1999), the scenario 
changes at z 3.5. At this redshift, the QSO contribution 
of ionizing photons falls below the critical limit needed 
to photoionize the IGM; and by z=5, it will fall short of 
the critical value by a factor of ^4. This implies that at 
high redshift, the contribution from young stars may be- 
come the dominant contributor to the background, with 
the caveat that the space density of star-forming galaxies 
would have to be maintained at the level observed at z « 
3, and that most of their UV photons would have to be 
free to escape into the IGM. The Devriendt et al. (1998) 
models lead to the conclusion, however, that the galaxy 



contribution to the UV background is negligible at high 
redshifts. 

In conclusion, the proximity effect data at present reflect 
that the UV background at 2 < z < 4 is QSO dominated. 
The discrepancies between this model at low and high red- 
shifts (Kulkarni & FaU 1993, Williger et al. 1994) indicate 
that the contribution to the background from galaxies may 
be of larger relative importance. We plan to undertake an 
analysis of the proximity effect at low redshifts from a large 
sample of QSO spectra taken with the Faint Object Spec- 
trograph on the Hubble Space Telescope to place better 
constraints on the background at 0.5 < z < 2. Further ob- 
servations of objects at z > 4 are also of particular interest 
to this subject. 
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Fig. 1.— 



Spectrophotometry of z f» 2 QSOs; Dashed Une indicates the power law continuum fit; dotted line indicates the Icr errors 
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Fig. 3. — Lyman limit luminosity versus redshift for the proximity effect dataset; squares- QSOs from which low redshift line sample was 
taken; crosses- QSOs from which high redshift line sample was taken; the line marks the boundary between low and high luminosity QSOs 
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Fig. 5. — x'^ of binned data with respect to the ionization model with a constant J{uo) versus log[J(i^o)]: (a) DB96 sample; (b) all lines, 
hya QSO redshifts; (c) all lines, using narrow line redshifts where available, Lya redshifts otherwise; (d) all lines, narrow line redshifts where 
available, hya redshifts + 400 km s^^ otherwise; (e) high z lines, QSO redshifts as in case (d); (f) low z lines, QSO redshifts as in case (d); 
(g) high luminosity QSOs, QSO redshifts as in case (d); (h) low luminosity QSOs, QSO redshifts as in case (d); (i) weak lines only; 0.16 
A <W < 0.32 A, QSO redshifts as in case (d); (j) lines from QSOs with damped Lya systems only, QSO redshifts as in case (d) 
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Fig. 6. — Number distribution per coevolving rcdshift coordinate for the best fit values of J(i'o) (BDO method); (a-j) same as Fig. 
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Fig. 8. — Number distribution per coevolving rcdshift coordinate for the best fit values of J(i'o) (KF method); (a-f) same as Fig. 
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Fig. 9. — Histograms of redshift differences with respect to the [OIII] A5007 Une (a) Lya; (b) Mg II; (c) Balmer lines ; dotted lines- data 
of Laor et al. (1995), dashed lines- data of Nishihara et al. (1997) 
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Fig. 10. — Sample simulation spectra plotted with data, flux scale is arbitrary. (a)0014+813: (i)data, (ii)input log[J(i'o)]=-23.0, (iii)input 
log[J(i^o)]=-22, (iv)input Iog[J(i/o)]=-21.3, (v)input log[J(i^o)]=-20.0, (vi)input log[J(i'o)]=-19.0; (b)1700+643: (i)-(vi) same as in (a) 
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Fig. 11. — Number distribution per coevolving redshift coiirdinate for the best fit values of J{uo) listed in Table [?[; solid lines- simulation, 
dotted lines- data, scaled by the relevant value of A.o in TableM: (a)input log[J(i/o)]=-23.0; (b)input log[ J(fo)]=-22; (c)input log[J(i^o)]=-21.3; 
(d)input log[J(i/o)]=-20.0; (e)input log[J(i^o)]=-19.0 
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Fig. 12. — Curve of growth effects: ratio of proximity effect to nonproximity effect rest equivalent widtlLfor all lines versus nonproximity 
effect rest equivalent width; solid line represents the detection threshold WpE= 0.32 A; (a)-(e) as in Fig. Ill 
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Fig. 13. — Measurements of log[J(!-'o)] vs. redshift: points and error bars are taken from Table |8[ The upper limit set by Bunker et al. 
(1998) at z~3 is included. Measurements over extended redshift ranges and the errors in those measurements are indicated by boxes. The 
solid curves are derived from the Haardt & Madau (1996) model for the HI photoionization rate as a function of redshift for QSO spectral 
indicies of (lower curve) and 2 (upper curve). Overall, measurements at z=2-3 agree well with one another and with the predictions of the 
Haardt & Madau (1996) model. 
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Fig. 14. — Power law fits to log[J(i/o)] as a function of redsiiift: J(!-'o,z) = J(fo,0)(l + z)^ . The dashed lines indicate the two lowest 
fits lathe data: (j,log[J(i/o, 0)])= (5.12,-23.97) and (-4.16,-18.76). The solid curves are the Haardt & Madau (1996) models as shown in 
Figure n3. The Haardt & Madau (1996) models are turning over at the redshift of the data, precluding a strong constraint on the parameters 
j and Jl^OjO); but the lowest fit extends to low redshift to match the Kulkarni & Fall (1993) measurement shown by the box, while the 
next lowest fit extends to high redshift to match the Williger et al. (1994) measurement, the point at z=4.5. 
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Table 1 

Spectrophotometry Observations 



QSO 


Dato 


Exposure 


Airiiiass 


Wa vol c 1 igt h C o vo r age 


0006+020 


29Novl994 


1800 


1.15 


3150-6385 


0027+018 


22Sepl992 


1800 


1.28 


3467-6475 


0037-018 


29Novl994 


2400 


1.27 


3150-6385 


0049+007 


29Novl994 


1800 


2.05 


3125-6380 


0123+257 


29Novl994 


1800 


1.55 


3125-6380 


0153+744 


29Novl994 


1800 


1.38 


3125-6380 


0348+061 


22Sepl992 


1800 


1.13 


3465-6475 


1323-107 


28Marl995 


1800 


1.56 


3115-6400 


1346-036 


28Marl995 


1800 


1.27 


3115-6400 


1422+231 


22Aprl996 


1800 


1.31 


5235-7554 


1603+383'' 


04Julyl995 


450 


1.03 


3663-7544 


2134+004 


22Sepl992 


1800 


1.29 


3465-6483 


2251+244 


29Novl994 


1800 


1.01 


3150-6385 


2254+022 


22Sepl992 


1800 


1.18 


3470-6480 



Instrument Set-up for 1422+231: SO B&C, 600 1 mm"! I''* order, A6=6681 A, 1.5" slit; for 1603+383: FLWO FAST, 300 1 
mm-l 1"* order, Ai,=4750, 3" slit 

^spectrum donated by Hamburg/CfA Bright Quasar Survey in advance of publication 



Table 2 

Summary of Narrow Emission Line Observations 



Name 


V 


Instrument 


Date 


Exposure (sec.) 


Wavelength Coverage (/im) 


0000-263 


17.5 


OSIRIS 


27Jull994 


4800 


1.20 - 


1.46 


0014+813 


16.5 


CRSP 


07Decl993 


1200 


1.18 - 


1.26 






FSpec 


26Novl993 


5280 


1.96 - 


2.39 


0114-089 


17.1 


CRSP 


04Decl993 


3180 


1.10 - 


1.35 


0636+680 


19.0 


CRSP 


04Decl993 


4800 


1.09 - 


1.35 






CRSP 


05Decl993 


2820 


1.96 - 


2.10 


0836+710 


16.5 


B&C 


29Marl995 


1800 


0.52 - 


0.96 


0956+122 


17.5 


CRSP 


04Decl993 


8220 


1.10 - 


1.35 






FSpec 


27Novl993 


5280 


1.96 - 


2.38 


1159+124 


17.5 


CRSP 


05Decl993 


3180 


1.09 - 


1.35 






FSpec 


29Novl993 


4320 


1.97 - 


2.38 


1207+399 


17.5 


FSpec 


21Mayl994 


600 


1.98 - 


2.41 


1208+101 


17.5 


CRSP 


06Decl993 


4800 


2.00 - 


2.42 


1408+009 


18.0 


FSpec 


02Aprl996 


3840 


1.46 - 


1.73 








02Aprl996 


1920 


1.99 - 


2.40 


1422+231 


16.5 


FSpec 


21Mayl994 


1920 


1.98 - 


2.41 


1435+638 


15.0 


FSpec 


02Aprl996 


1920 


1.99 - 


2.40 


2126-158 


17.3 


CRSP 


05Decl993 


3180 


1.08 - 


1.35 






OSIRIS 


24Sepl994 


7680 


1.96 - 


2.35 
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Table 3 

Spectrophotometric Properties 



QSO 


Nhi (lO^o cm-2) 


£^'"'(912 A) 


a 


tobs 


f^(912 A) 


Ref. 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


0006+020* 


3.02 




0.26 


354 (1450 A) 


313 


1 


0027+014 


2.93 




-0.38 


219 (1450 A) 


183 


1 


0037-018 


2.81 




-0.27 


45 (1450 A) 


51 


1 


0049+007 


2.67 




0.31 


324 (1450 A) 


280 


1 


0123+257* 


6.88 




1.12 


237 (1450 A) 


141 


1 


0150-202* 


1.29 






529 (1430 A) 


430 


2,3 


0153+744* 


22.74 




0.18 


1023 (1450 A) 


940 


1 


0226-038 


2.35 






582 (1800 A) 


425 


4 


0348+061 


12.33 




0.12 


513 (1450 A) 


485 


1 


0400+258 


7.82 




1.54 


199 


5 


0747+610 


4.77 






500 (1800 A) 


365 


4 


0819-032 


6.16 




0.33 


63 (1450 A) 


54 


6 


0836+710* 


2.93 








652 




0848+155 


3.14 




0.07 


198 (1450 A) 


191 


7,8 


0936+368 


1.36 








386 




0952+335 


1.37 








370 




0955+472* 


1.04 








188 




0956+122 


3.10 


140 


0.49 


448 (1450 A) 


356 


9 


1009+299 


2.30 








1217 




1207+399 


2.10 




0.59 


319 (1450 A) 


242 


1,8 


1210+175* 


2.67 






285 




1231+294 


1.54 








980 




1323-107 


2.64 




-0.30 


303 (1450 A) 


349 


1 


1329+412* 


0.99 




0.33 


750 


10 


1337+285* 


1.17 








339 




1346-036 


2.51 




0.091 


458 (1450 A) 


439 


1 


1358+115* 


1.81 




1.10 


345 (1450 A) 


207 


6 


1406+492 


1.77 








392 




1408+009 


3.04 




0.91 


99 (1450 A) 


64 


1 


1421+330 


1.23 


58 


0.54 


914 (1450 A) 


711 


11,7 


1422+231* 


2.52 




-1.21 


211 (1450 A) 


371 


1 


1435+638 


1.68 


55 




1244 (1800 A) 


909 


12,4 


1603+383* 


1.32 




0.36 


550 (1450 A) 


464 


1,13 


1604+290 


3.24 






428 




1715+535 


2.69 


36 


1.26 


875 (1800 A) 


371 


11,10,4 


2134+004 


4.03 




0.04 


35 (1450 A) 


34 


1,14 


2251+244* 


5.18 




1.53 


243 (1450 A) 


119 


1 


2254+024 


5.32 




0.20 


116 (1450 A) 


106 


1 


2310+385 


10.62 






419 




2320+079 


5.04 








306 




2329-020 


4.45 








451 





^QSO name; an asterisk denotes a metal line system within 5000 km s ^ of the QSO emission redshift; in the case of 1422+231, the 
QSO is a known lens 

'^Galactic ^hi in units of 10^*^ cm^^ from program GOLDEN using Stark et al. (1992) 
'^Observed flux in /ijy at tlic Lyman limit from reference in (g) 

'^Observed spectral index between Lya and C IV emission lines or in the vicinity of the flux listed in (e) from reference in (g); in 
general, values are based upon spectra corrected for Gala<;tic reddening if E(B-V) > 0.03 

^Observed flux in /njy at the rest wavelength indicated in parentheses from reference in (g) 

Extrapolated Lyman limit flux in ^Jy from measured flux in (e), when available, or V magnitude given in Table 1 of Paper I.; if no 
observed spectral index available, value of 0.46 used (Francis 1996) 

SrEFERENGES: (l)this paper; (2)MacAlpine & Feldman 1982; (3)GrifHth et al. 1994; (4) Steidel & Sargent 1991; (5)Gheng,GaskeU, 
& Koratkar 1991; (6)Pei, Fall & Bechtold 1991; (7)Uomoto 1984; (8)Barthel et al. 1988; (9) Sargent, Steidel, & Boksenberg 1989; 
(lO)Baldwin, Wampler, & GaskeU 1989; (ll)Koratkar, Kinney, & Bohlin 1992; (12)Lanzetta, Turnshek & Sandoval 1993 (13) Ham- 
burg/GfA Bright Quasar Survey (unpublished); (14)Perez, Periston, & Moles 1989 



37 



Table 4 
Measurements of J(i'o) 



Sample 
(a) 


lines 

(b) 


7, norm. 

(c) 


method 
(d) 


log[(J(i^o)] 


(f) 


■^points 

(g) 


(h) 


Figure 
(i) 


1 


518 


1.9260,5.8882 


BDO 


-21-40tJ-ig 


3.05 


7 


0.88 


6(a) 


1 


518 


1.9260,3.9709 


ML 


-21.58«:^,» 


20.3 


6 


0.0024 


8(a) 


2 


1286 


1.6749,7.5723 


BDO 


-20.90«;« 


5.22 


7 


0.63 


6(b) 


2 


1286 


1.6749,4.6637 


ML 


-20.83tojg 


6.32 


6 


0.38 


8(b) 


3 


1286 


1.6749,7.5723 


BDO 


-21 no"'"'''^^ 

Z1.UU_Q 3g 


7.19 


7 


0.40 


6(c) 


3 


1286 


1.6749,4.6709 


ML 


-20.83«:^* 


7.41 


6 


0.28 


8(c) 


4 


1286 


1.6749,7.5723 


BDO 


21 15+0-17 


6.54 


7 


0.47 


6(d) 


4 


1286 


1.6749,4.6617 


ML 


-21 17+019 


3.53 


6 


0.73 


8(d) 


5 


763 


-0.2848,110.13 


BDO 


on 7c;+0-l6 


3.31 


7 


0.85 


6(e) 


5 
6 


763 
523 


-0.2848,69.934 
1.3754,10.240 


ML 
BDO 


-21.18t°i? 
21 15+011 


4.92 
3.97 


5 
7 


0.42 
0.78 


8(e) 
6(f) 


6 


523 


1.3754,7.4759 


ML 


-21 46+0-3" 
Zl.w_|3 29 


15.5 


6 


0.016 


8(f) 


7 


261 


2.3284,2.6809 


BDO 


-21 45+0-« 


3.03 


7 


0.88 


6(j) 


8 


666 


1.5361,9.1237 


BDO 


01 ocr-t-0.28 
-zl.zo_Q 45 


4.32 


7 


0.74 


6(g) 


9 


620 


2.0242,4.6980 


BDO 


-Zi.UO_Q 42 


3.49 


7 


0.83 


6(h) 


10 


671 


0.5468,24.655 


BDO 


-20 4.5+0-37 


3.05 


7 


0.88 


6(i) 



^(1) DB96 sample; (2) all lines, Lya QSO redshifts; (3) all lines, using narrow line redshifts where available, 
hya redshifts otherwise; (4) all lines, narrow line redshifts where available, Ly« redshifts -|- 400 km s~l 
otherwise; (5) high z lines, QSO redshifts as in case (4) ; (6) low z lines, QSO redshifts as in case (4); (7) 
lines from QSOs with damped Lyo systems only, QSO redshifts as in case (4); (8) high luminosity QSOs, 
QSO redshifts as in case (4); (9) low luminosity QSOs, QSO redshifts as in case (4); (10) weak lines only: 
O.ieA <W < O.32A, QSO redshifts as in case (4) 

''number of Lyo forest lines in sample 

"^Equ. 1^ parameters 7 and Ao from maximum likelihood fit to data; when the method listed is ML, the 
normalization listed is equal to AoiNii^/No)^~^{f} — l)~i (see text, Paper I) 

'^BDO- BajtUk, Duncan, & Ostrikcr (1988), ML- maximum likelihood, see Kulkarni & Fall (1993) 
®Best fit value of log[{J (uq)] in units of ergs s~i cm"^ Hz~i sr^i 

of data vs. the ionization model used 
^number of points used to calculated 

probability for the ionization model used 
'Figure displaying number distribution per coevolving redshift interval, dj\f /dX-y 
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Table 5 

Emission Line Redshifts for QSOs Included in Measurement of 



QSO 


z 


line'' 




0000-263* 


4.111 


Lya 


1 




4.116 


Mg II 


2 


0001+087 


3.243 


Lya 


3 


0002+051 


1.899 


Ly« 


4 




1.899 


Mg 11 


5 


0002-422 


2.763 


Lya 


6 


0006+020* 


2.340 


Lya 


2 


0014+813 


3.386 


Lya 


1 




3.379 


Mg II 


2 




3.404 


H/3 


2 


0027+014 


2.333 


Lya 


2 




2.310 


H/3 


7 


0029+073 


3.261 


Lya 


1 


0037-018 


2.341 


Lya 


2 


0049+007 


2.275 


Lya 


2 




2.279 


[OIII] A5007 


8 


0058+019 


1.959 


Lya 


9 




1.964 


Mg II 


5 


0100+130 


2.690 


Lya 


6 


0114^089 


3.194 


Lyo 


9 




3.192 


Mg II 


2 


0119-046* 


1.951 


Lya 


1 




1.964 


Mg II 


5 


0123+257* 


2.358 


Lyo 


10 




2.370 


[OIII] AA4959,5007 


8 


0142-100* 


2.727 


Lya 


9 


0150-202* 


2.148 


Lya 


2 




2.149 


Mg II 


5 


0153+744* 


2.340 


Lya 


2 




2.341 


[OIII] A5007 


8 


0226-038 


2.067 


Lya 


2 




2.073 


Mg II 


5 




2.073 


[OIII] AA4959,5007 


8 


0237-233 


2.224 


Lya 


9 




2.200 


[OIII] A5007 


7 


0256-000 


3.374 


Lya 


1 


0301-005 


3.228 


Lya 


1 


0302-003* 


3.286 


Lya 


1 


0334^204 


3.131 


Lya 


1 


0348+061 


2.057 


Lya 


2 




2.056 


Mg II 


5 


0400+258 


2.108 


Lya 


2 


0421+019 


2.050 


Lya 


4 




2.056 


Mg II 


5 


0424^131* 


2.165 


Lya 


11 




2.166 


Mg II 


5 




2.163 


Ha 


12 


0453-423 


2.656 


Lya 


6 


0636+680 


3.167 


Lya 


1 




3.184 


Mg II 


2 




3.187 


H/3 


2 


0731+653 


3.033 


Lya 


1 


0747+610 


2.491 


Lya 


2 


0831+128 


2.739 


Lya 


3 


0836+710* 


2.189 


Lya 


2 




2.197 


Mg II 


2 



Table Tj Coutiuucd 



QSO 


z 




line"^ 






9 918 
Z.Zlo 


rr^TTTl \c;nn7 

[(^ilij AuUUi 




nQQVj-.! no 


Q Q9Q 


Lya 


Q 

y 




9 ni Q 
z.uiy 


Lya 


9 

z 




o ni y1 
Z.U14 


A^fT TT 
IVlg 11 


o 


HQ/I AQ* 


1 Q9C^ 

i.yzo 


Lya 


lo 




1 Q99 

j..yzz 


\l[rT TT 
IVlg 11 


O 


uyuo-rioi 


^ 1 7*^ 


Lya 


O 


uyio-rU t z 


9 78^; 
Z. 1 oO 


Lya 


Q 

y 


uyoD+oDcS 


o nor; 
z.UzO 


Lya 


z 


uyoo-riiy 


1 Q9 
o. lyz 


Lya 


O 




9 K(\A 
Z.OU4 


Lya 


z 


uyoo-r^ J z 


9 4S9 


Lya 


9 
Z 


flQ'^fiJ-l 99 

uyoD-r A.^-^ 


o.Uoo 


Lya 


9 
Z 




9QQ 

o. zyy 


A/Trr TT 
iVlg 11 


9 
Z 




Q Q1 A 




z 




o.oUo 


rr^TTTl \'=;nn7 

[t^lllj AOUUi 


z 




z.uoo 


Lya 


2 


iUi / -\-z.o\) 


1 Q9S 

1 .yzo 


Lya 


Q 

y 


iUoo-|-lo ( 


no9 
o.uyz 


Lya 


o 


111 c;_i_nfin* 


1 797 
i. ( Z i 


Lya 


1 /I 

14 


1 1 f^Qj-l 94* 


O.OUO 


Lya 


Q 

y 




o.OUo 


AAfT TT 
iVlg 11 


9 
Z 




o.4y / 


Hp 


Z 




407 

o.4y ( 


rOTTTl \ ^f\f\7 
[WlllJ AOUU ( 


9 
Z 


1 9nfi_i_i 1 o 
izUD-|-iiy 


o.iUo 


Lya 


Q 

O 


izu 1 "Toyy 


9 /I C;i 
Z.40i 


Lya 


9 
Z 






TT/-W 
llCx 


9 
Z 




fi99 
vS.oZZ 


Lya 


Q 
O 




o.ooo 


TT R 


9 
Z 




Q Qn9 


rr^TTTl \c;nn7 

[^^lilj AOUUi 


9 
Z 


1 91 n_i_i 71^* 


Z.004 


Lya 


9 
Z 




Z.DUD 


Lya 


1 4 
11 


1 99c^ m 1 


9 QQ1 
Z.cSoi 


Lya 


10 


1 99c;_L'i1 7 


2 200 


Lya 


o 




Z.zzO 


rrMTTl \ c;nn7 
[UlilJ AoUU( 




1 9Q1 _L.9Q/1 

iz.3i-T-zy4t 


9 ni Q 


Lya 


9 
Z 




9 flAI 
Z.U41 


Lya 


Q 

y 


1 Q 1 _l_ y1 70 


o Kon 

z.oyu 


Lya 


Q 

o 


1323-107 


2.360 


Lya/C IV 


16 


1329+412* 


1.934 


Lya 


2 


1334^005 


2.842 




3 


1337+285* 


2.541 


Lya 


2 


1346-036 


2.356 


Lya 


2 




2.368 


Mg II 


12 




2.362 


[OIII] A5007 


8 




2.367 


Ha 


12 


1358+115* 


2.589 


Lya 


2 


1400+114 


3.177 


Lya 


3 


1402+044 


3.208 


Lya 


17 


1406+492 


2.161 


Lya 


2 


1408+009 


2.262 


Lya 


2 




2.260 


[OIII] A5007 


2 




2.265 


Ha 


2 


1410+096 


3.313 


Lya 


3 


1421+330 


1.903 


Lya 


2 




1.906 


Mg II 


5 
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Table 5 — Continued 



QSO 


z 


line'" 


Ref.'= 


1422+231* 


3.624 


Lya 


2 




3.630 


H/3 


2 




3.623 


[OIIIJ A5007 


2 


1435+638 


2.063 


Lyo 


2 




2.061 


Mg II 


5 




2.066 


[OIII] AA4959,5007 


8 




2.065 


Ha 


2 


1442+101 


3.560 


Lya 


17 


1451+123 


3.251 


Lya 


3 


1511+091* 


2.877 


C IV 


9 


1512+132 


3.120 


Lya 


3 


1548+092 


2.759 


Lya 


9 


1601+182* 


3.227 


Lya 


3 


1602+178* 


2.989 


Lya 


3 


1603+383* 


2.510 


Lya 


18 


1604+290 


1.962 


Lya 


2 


1607+183 


3.120 


Lya 


17 


1614+051 


3.216 


Lya 


17 




3.214 


[OIIIJ A5007 


19 


1623+269* 


2.526 


Lya 


3 


1700+642 


2.744 


Lya 


15 


1715+535 


1.935 


Lya 


2 




1.932 


Mg II 


5 


1738+350 


3.239 


Lya 


3 


1946+770 


3.020 


Lya 


15 


2126-158 


3.280 


Lya 


6 




3.284 


Mg II 


2 




3.298 


H/3 


2 




3.292 


[OIII] A5007 


2 


2134+004 


1.941 


Lya 


2 


2233+131 


3.301 


Lya 


1 


2233+136 


3.207 


Lya 


1 


2251+244* 


2.335 


Lya 


17 




2.359 


[OIII] A5007 


8 


2254+024 


2.089 


Lya 


2 




2.090 


Mg II 


5 


2310+385 


2.179 


Lya 


2 




2.181 


[OIII] AA4959,5007 


8 


2311-036 


3.041 


Lya 


1 


2320+079 


2.088 


Lya 


2 


2329-020 


1.896 


Lya 


2 



''Objects witli both [OIII] and Lya, Mg II, or Balmer redshifts were used 
to construct histograms in Figure ^; Objects marked with an asterisk are ex- 
cluded from the proximity effect analysis on the basis of associated absorption 
or gravitational lensing 

'^Emission lines used to measure redshift. 

REFERENCES: (1) Sargent et al. 1989; (2) this paper; (3) B94; (4) Young et al. 1982a; 
(5) Steidel & Sargent 1991; (6) Sargent et al. 1980; (7) Baker et al. 1994; (8)M'^Intosh et al. 
1999a (**1225+317 measurement quoted as uncertain due to weak [O III] emission and low 
S/N); (9) Sargent et al. 1988; (10) Schmidt & Olsen 1968; (11) Burbidge 1970; (12) Espey 
et al. 1989; (13) Young et. al 1982b; (14) Wills & Wills 1979; (15) DB96; (16) Kunth et al. 
1981; (17) Barthel et al. 1990; (18) Hamburg/CfA Bright Quasar Survey (unpubhshed); 
(19) Bremer & Johnstone 1995 
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Table 6 
Ionization Rates 



A,B,Zc,S 


Ref. 




%2 




(a) 


(b) 


(c) 


6.7e-13 s-1,0.43,2. 30,1.95 


1 


11.2 


0.12 


5.6e-13 s-1,0.60,2.22,1.90 


2 


11.8 


0.10 


1.2e-12 s-1,0.58,2. 77,2.38 


2 


7.15 


0.41 



'*(l)Haardt & Madau (1996); (2)Fardal et al. (1998) 

of data vs. the BDO ionization model 
'^X^ probability for the BDO ionization model 



Table 7 
Simulation Results 



Input log[{J{uo)] 
(a) 


(b) 


log[(J(!/o)] recovered 
(c) 


(d) 




-23.0 
-22.0 
-21.3 
-20.0 
-19.0 


1.5722,11.043 
1.6869,8.8367 
2.6267,2.6960 
2.2511,3.8084 
2.0302,5.2704 


OO 7c:+0.28 
"^^•'^-0.19 

-21.80«:« 
-21 OO"*""'^* 

-18.50+';-J6 


11.2 
11.3 
2.68 
3.90 
5.28 


0.12 
0.12 
0.91 
0.79 
0.62 



^value of log[( J(i^o)] used for modifying absorber column densities according to Equations^ and 

^Equ. 1^ parameters 7 and ^0 from maximum lilselihood fit to data 

'^value of log[(J(j'o)] from simulated spectra using the standard BDO technique 

'^X^ of data vs. the BDO ionization model 

probability for the BDO ionization model 



Table 8 

Literature Proximity Effect Measurements of J(i'o) 



log[(./(!-0)] 


z 






Ref. 


on 9+0.8 

^■"■■^-o.e 


0.16-0. 


,99 


13 


1 


-20.5+f'3 


1.8-2. 


,3 


3 


2b 


-21.1±d.6 


2.0 




1 


2 


-21.15 


3.2 




1 


3 


-21.3 


3.6 




1 


4 


-21.0+0.5 


1.7-3. 


,8 


38 


5 


-21.0+0.5 


1.7-3. 


,8 


19 


6 


-20.5 


1.6-4. 


,1 


49 


7 


91 0+0. 08 


1.7-4. 


1 


10 


8 


91 1+0.15 

-^^■'^-0.27 


1.7-4. 


1 


74 


9 


-21 0+°-" 
^'-■"-0.16 


2.0-4. 


,5 


11 


10 


-22.0 - -21.5 


4.5 




1 


11 



* (1) Kulkarni & Fall 1993; (2) Fernandez-Soto et al. 1995; (3) Giallongo et al. 1993; (4) Cristiani et al. 1995; (5) BDO: (6) 
LWT; (7) B94; (8) Giallongo et al. 1996; (9) this paper; (10) Cooke et al. 1997; (11) Williger et al. (1994) 



'measured from the proximity effect to due a foreground QSO; not able to set an upper limit 



